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ABSTRACT 

Calculations carried out to model the evolution of Star 32 under different assumptions 
about the stellar wind mass-loss rate provide robust limits on the present mass of the 
star. The obtained range is 31 to 35.5 M©, which is in very good agreement with the 
orbital solution of Orosz et al., namely 28.3 to 35.3 Mq. The initial mass of Star 32 
had to be in the range 35 to 40 and the present age of the system is 3.7 to 4.0 
Myr. 
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1 INTRODUCTION 

LMC X-1 was one of the first X-ray binaries discovered, 
and the first one found in Magellanic Clouds (Mark et al. 
1969). The X-ray source is persistent, although variable, 
the luminosity varying by one order of magnitude (Liu, van 
Paradijs & van den Heuvel 2005). The source is very bright 
(Lx ~ 1-^2x10^* erg s"^ Johnston, Bradt & Doxsey 1979). 
Its X-ray spectrum was found to be relatively soft (fcT ~ 2.7 
keV, Markert & Clark 1975), which prompted Hutchings, 
Crampton & Cowley (1983, hereafter H83) to notice that 
it was similar to the two then known black hole candidates: 
Cyg X-1 and LMC X-3. The optical identification was firmly 
determined only quite recently (Cowley et al. 1995). The dif- 
ficulty with the identity of the optical counterpart was that 
there were two comparably good candidates. Initially the 
bright (V — 12.0) B5 supergiant R148 was favored as the 
counterpart (Jones, Chetin & Liller 1974; Rapley & Tuohy 
1974; Johnston et al. 1979). The other, less favored candi- 
date was the fainter {V ~ 14.5) star denoted by Cowley et 
al. (1978) as star no. 32. This star was observed spectro- 
scopically by PakuU (1980), who noticed that the spectra 
"show He II 468.6- nm and N III-C III 464-465- nm emissions 
with strength comparable to that seen in most massive X- 
ray binary systems" . Both optical candidates were observed 
spectroscopically by Hutchings et al. (1983). They found no 
detectable periodic variability for R148, but radial velocity 
measurements of Star 32 have shown it to be a binary with 
an orbital period of approximately 4 days. Two equally good 
fits to the observations were obtained for values of the pe- 
riod equal to 3.908 and 4.038 d. The authors also measured 
the velocities of emission lines, discovered by PakuU, which 
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were found to vary in antiphase with the lines of Star 32. 
These lines probably originated near the secondary (com- 
pact) component, and, from the relative amplitudes, the au- 
thors deduced that the mass ratio (defined as the mass of 
the optical to the mass of the compact component) should 
exceed 2. The estimated most likely masses of the compo- 
nents were about ~ 14 and ~ 4 Mq, which suggested that 
the secondary may be a black hole. White and Marshall 
(1984) noticed that the very soft X-ray spectrum of LMC 
X-1 places it close to other black hole candidates such as 
LMC X-3 and Cyg X-1 in an X-ray colour-colour diagram. 
After further optical spectroscopy, Hutchings et al. (1987, 
hereafter H87) refined the orbital period to 4.2288 d, and 
the masses of the components to 20 and 6 M0. The evi- 
dence that Star 32 is indeed the optical counterpart of LMC 
X-1 looked quite strong, but it was firmly established only 
after Cowley et al. (1995) compared the position of Star 32 
with the position of LMC X-1 from ROSAT observations. In 
2006, Levine & Corbet found the X-ray orbital period from 
data from RXTE. Their period was equal to 3.9081 ± 0.0015 
d, which was somewhat shorter than the then accepted H87 
period (4.2288 d), but was equal to one of the original (1983) 
propositions of H83. 

In 2005 I calculated the preliminary evolutionary tracks 
modelling the evolution of Star 32 (Ziolkowski 2006). I found 
that the present mass of this star should be in the range 
24 to 33 Mq. However, this result was difficult to reconcile 
with the then available spectroscopic data. The reason was 
as follows. 

We can estimate the radius iiopt of the optical com- 
ponent from the observationally determined luminosity and 
the effective temperature. For spectral type 07 HI (H83) 
we have Tc ^ 36000 K and B.C. « -3.50. H83 estimated 
reddening as Eb-v ~ 0.35, which implies Av « 1.1. With 
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V « 14.8 and the distance modulus to the Large Magel- 
lanic Cloud (LMC) equal to 18.5 we have Afboi ~ —8.3 or 
log(L/Lo) « 5.2. With this luminosity and effective tem- 
perature, one obtains i?opt ~ 10.5i?o- 

If we assume the mass of the optical component Mopt 
to be a free parameter, then after selecting its value we can 
use two equations to solve for the inclination of the orbit i 
and the mass ratio q = Mopt/Mx (where Mx is the mass of 
the compact component). One of these equations makes use 
of the mass function: 

/(Mx) = Moptsin^VkCl + Qf]- (1) 

The other relates the radius of the star to the size of the 
orbit: 

Ropt = Rri^ X /rl = 

= /RL(0.38-|-0.21og g)A = 

= /RL(0.38 + 0.21og g)ai(l + g), (2) 

where -Rrl is the radius of the Roche lobe (e.g. Paczyriski 
1971) around Star 32, /rl is the fill-out factor (/rl = 
Ropt/RnL), A is the orbital separation of the binary com- 
ponents and oi is the radius of the orbit of Star 32. 

Adopting, after H87, the (then accepted) values of the 
orbital period (4.2288 d) and of the radial velocities am- 
plitude {Kopt = 69 km/s), one has /(Mx) = O.144M0 and 
aisin i = 5.76R0 for the mass function and the radius of the 
orbit of the optical component, respectively. 

Inserting these data, equations (1)— (2) can be written 

as 

Moptsm^i/[q{l + qf] = 0.144, (3) 

Ropt = /rl(0.38 + 0.2 log q){l + q) x 5.76/sin i. (4) 

Recall that both H83 and H87 indicate that the surface 
of the star must be near the Roche limiting surface (/rl ^ 
0.9). Then, solving equations (3) — (4) for i and q, it is easy 
to verify that the observational condition q > 2 (H83) can 
be satisfied only for Mopt < 8 M© if /rl = 0.95 or for 
Mopt ^ 9.5 M0 if /rl ~ 0.9. Allowing for observational 
errors, one might increase these upper limits slightly for the 
value of Mopt, but this would not bridge the gap between 
the orbital solution (Mopt < 8 9 Mq , as shown above) and 
the evolutionary constraints (Mopt > 24 M©). 

Therefore, I concluded my 2006 considerations with the 
statement "There are only two possible solutions of this 
discrepancy: either observations of Star 32 (spectroscopy 
and/or photometry) arc in serious error or Star 32 is es- 
sentially a helium star (with only a small remnant of the 
hydrogen rich cnvcloijc)'" . 1 concluded that future observar 
tions should solve this problem. 



2 NEW ORBITAL SOLUTION 

Later observations indeed solved the problem. Orosz et al. 
(2009) published a detailed analysis of their new extensive 
spectroscopic and photometric observations of Star 32 to- 
gether with a refined analysis of the ASM data from RXTE 
observations of LMC X-1. Their new orbital solution greatly 
improved the earlier rough estimates, removed the puzzles 
(mass versus luminosity of the optical component) and pro- 
duced a fairly consistent picture of the binary system. First, 



they found that the optical orbital period is not 4.2288 d 
(as was generally accepted after H87), but rather 3.90917 
± 0.00005 d, in agreement with the X-ray period suggested 
by Levine & Corbet (2006), and also with one of the two 
original orbital period candidates of H83. After a thorough 
analysis, they determined the masses of the components as 
Mopt = 30.62 ± 3.22 M® and Mx = 10.30 ± 1.34 M©. They 
also determined quite precisely the luminosity and the ef- 
fective temperature of Star 32 as log(L/Lo) = 5.50 ± 0.05 
and Tc = 33200 ± 500 K. They found that these values are 
consistent with the core burning star of initial mass ^ 35 

Mq. 

It thus became clear that the reasons for the dramatic 
discrepancy discussed in Ziolkowski (2006) were: (i) a serious 
(by a factor of 2) underestimate of the luminosity of Star 32 
and (ii) an overestimate (by ~ 8.5 per cent) of its effective 
temperature. These two factors led to a serious (by ~ 40 
per cent) underestimate of its radius (~ 10.5 R©, instead 
of ^ 17 R0). This small value for the radius of Star 32, 
together with the requirement that it should nearly fill its 
Roche lobe, led to an uncomfortably small mass of this star. 

Note that even the underestimated luminosity of 
log(L/L0) ~ 5.2 was much too high for the mass of Star 32 
proposed by 1183 (~ 14 M©; and this relatively high value of 
the mass was achieved only by stretching the then observed 
value of radius from ~ 10 Rq to ~ 12 Rq). 

After new precise data concerning the LMC X-l/Star 
32 binary system became available, I decided to carry out 
new evolutionary calculations to determine more precisely 
the past evolution and the present evolutionary state of Star 
32. 



3 EVOLUTIONARY CALCULATIONS FOR 
STAR 32 

3.1 General description 

I computed evolutionary tracks for core hydrogen-burning 
phase of stars with initial masses in the range 30-50 Mq . The 

Warsaw evolutionary code developed by Bohdan Paczjaiski 
and Maciek Kozlowski and updated by Ryszard Sienkiewicz 
and Alosha Pamyatnykh was used. An initial chemical com- 
position of X=0.7 and .Z=0.008, appropriate for LMC, was 
adopted. Opacity tables incorporating OPAL opacities (Igle- 
sias & Rogers 1996) as well as molecular and grain opacities 
(Alexander & Ferguson 1994) were used. The nuclear reac- 
tion rates are those of Bahcall & Pinsonneault (1995). The 
equation of state used was that of Livermore Laboratory 
OPAL (Rogers, Svenson & Iglesias 1996). Scmiconvcctivc 
mixing was neglected, as it is not important during the evo- 
lutionary phase considered (most of the models of interest 
had central hydrogen content Xc > 0.2). Similarly, any over- 
shooting at the border of the convective core was neglected, 
as this is even less important. 

The calculations were carried out under the assumption 
that the evolution starts from homogeneous configurations. 
This means that the consequences of the fact that some 
of the matter of the star, possibly dumped from the pro- 
genitor of the present black hole, could have modified the 
chemical composition, were neglected. It also means that 
the consequences of the fact that some nuclear evolution 



(hydrogen burning) could, possibly, have taken place while 
the mass of the star was smaller (prior to the mass transfer) 
were neglected as well. It seems that neither simplification 
significantly alters the outcome of the evolutionary calcula- 
tions. The orbital period is so short that any substantial 
mass transfer during past evolution seems unlikely, as it 
would have lead to a common envelope configuration and 
the merger of the two stars. Even if (which seems unlikely) 
there was significant mass transfer in the system, then it just 
reset the evolutionary clock and we can consider the evolu- 
tion of the mass gainer as starting anew from the zero-age 
main sequence (ZAMS), as a single star but with a higher 
mass. 
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Table 1. Parameters of selected evolutionary models for 
Star 32 



No. Mo Mopt /now logL -M 

[Mq] [Mq] [Lq] [lO-6M0/yr] 



1 


33.5 


32.34 


1 


5.446 


0.57 


2 


37.5 


35.93 


1 


5.545 


0.84 


3 


41.5 


29.66 


8.5 


5.451 


4.85 


4 


44.5 


33.66 


6 


5.551 


5.10 


5 


36 


31.28 


8.5 


5.443 


4.66 


6 


40 


35.19 


6 


5.550 


5.07 


7 


35 


31.79 


8.5 


5.446 


4.63 


8 


39 


35.58 


6 


5.550 


5.04 



3.2 Stellar wind mass loss 

The most uncertain element of the calculations of the early 
evolution of massive stars is the mass loss caused by the 
stellar wind. The uncertainty in the estimate of its rate is 
the single most important factor influencing the outcome 
of the calculations (see e.g. Ziolkowski 2005). The observa- 
tions seem to indicate that there is a substantial scatter of 
mass-loss rate among stars of similar luminosities and ef- 
fective temperatures. The commonly used formula derived 
by Hurley, Pols & Tout (2000, hereafter HPT), based on 
paramctrization of Nieuwenhuijzcn & de Jager (1990), gives 
the estimate of the mass-loss rate with an accuracy that is 
probably not better than a factor of two. Bearing this in 
mind, I introduced the multiplicative factor /sw applied to 
the HPT formula. In this way, the uncertainty in the theory 
of evolution could be, in some way, taken into account. 

The value of the factor /sw, which should be used to 
model successfully the evolution of a given star, might be 
quite different for seemingly similar stars. For some cases 
(e.g. HDE 226868, which is the companion to Cyg X-1), this 
value is close to 1 (Ziolkowski 2005). However, in some high 
mass X-ray binaries, the optical supergiant components are 
significantly under-massive for their luminosities (Ziolkowski 
1977). In some systems, such as Cen X-3, this undermassiv- 
ness is very serious and requires much stronger mass loss 
than the normal HPT stellar wind (Ziolkowski 1978). It 
seems that Star 32 is probably (as we shall see) closer to 
the case Cen X-3 than to HDE 226868, in this respect. 

The evolutionary calculations carried out to produce 
models reproducing the present evolutionary state of Star 
32 and using the HPT formula are quite successful (see Fig. 
la), except for one aspect: the value of the stellar wind mass- 
loss rates predicted by these models are too small by almost 
an order of magnitude when compared with the observations 
(see Table 1). The observed value M ^ -5 x 10"^ Mo/yr 
was determined by Orosz et al. (2009) from analysis of the 
orbital X-ray light curve. They modelled the light curve suc- 
cessfully, assuming that variability is due to the scattering of 
X-ray photons by electrons in the stellar wind over the vari- 
able (with the orbital phase) optical depth in the wind. This 
estimate agrees roughly with the one derived from the ob- 
served X-ray luminosity. It should be noted, however, that 
this estimate, being model dependent, is not as robust as 
the determinations of the parameters of the system and its 
uncertainty might be quite large. 

The observed value of M seems to be substantially (by 
a factor of ~ 6 - 8) larger than the value predicted by the 



NOTES: 

(i) No. denotes the number of the evolutionary sequence; 
sequences 1 and 2 were obtained using the HPT formula 
(see text); sequences 3 and 4 used formula (5); sequences 
5 and 6 used formula (6); and sequences 7 and 8 used for- 
mula (7) 

(ii) Mo denotes the initial (zero-age main sequence) mass 
of the optical component; Mopt, the mass at the evolution- 
ary phase corresponding to the present state of Star 32; 
/now , the multiplicative factor applied in formulae (5)— (7) ; 
other symbols have their usual meanings 



HPT formula for the present parameters of Star 32. It is 
clear, therefore, that to model the evolution of Star 32 (in- 
cluding its present mass-loss rate) successfully it is neces- 
sary to modify the HPT formula, increasing significantly 
the strength of the stellar wind; that is, to use values of the 
factor /sw that are substantially larger than 1. I have tried 
three such modifications: 

/sw = /now, (5) 
/sw = 1-1- (/now — ^){R — RzAMs)/{Rci — /?ZAMs), (6) 
/sw = 1 + (/now - l)[(i? - i?ZAMs)/(i?cr - i^ZAMs)]^ (7) 

In these equations, /now denotes the value of /sw cor- 
responding to the present state of Star 32, that is, the value 
one has to apply to reconcile the calculated mass flux of the 
stellar wind (for the evolutionary model corresponding to 
the present state of Star 32) with the observed one. It was 
found to be ~ 6 for the upper edge and ~ 8.5 for the lower 
edge of the observational error box (the crosses in Figs la- 
d). ii is the temporary stellar radius, J?zams is the initial 
value of the stellar radius (on the ZAMS) and Ra is the 
present radius of the inner Roche lobe around the optical 
component. 

The modification given by equation (5) is the simplest 
one, but is not very realistic. There is no reason why the 
stellar wind should be, from the very beginning, stronger 
than the typical one (HPT) by a factor as high as 6-8. 
However, there arc indications, dating from a long time 
ago (Ziolkowski 1977, Hutchings 1979), that the strength of 
the stellar wind from a component of a close binary system 
might increase significantly as the stellar surface approaches 
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Figure 1. The evolutionary tracks in the Hertzsprung- Russell diagram (solid lines). The dotted cross shows the observed position of 
Star 32. The number at the start of each tra<;k indicates the initial mass (in solar units). The numbers near the cross indicate the mass 
at the evolutionary phase corresponding to the present state of Star 32. Part (a) shows the tra<;ks obtained using the HPT formula to 
describe the stellar wind (see text); (b) shows tracks obtained using formula (5); (c) shows tra<;ks obtained using formula (6); and (d) 
shows tracks obtained using formula (7). The value of the parameter /now (see text) used for eax;h trax;k is given in Table 1. 



that of the Roche lobe. There is no quantitative descrip- 
tion of this effect available. The formuleie (6) and (7) try 

to incorporate it, in a crude way, by a linear or quadratic 
dependence on the distance from the Roche lobe. 

3.3 Results 

The evolutionary tracks in the Hertzsprung-Russell (HR) di- 
agram, obtained with the unmodified HPT formula and with 
formulae (5)-(7), are shown in Figs l(a)-(d). All tracks re- 
produce the present luminosity and effective temperature of 
Star 32 well, and the tracks obtained with formulae (5)-(7) 
reproduce, in addition, the present mass flux of the stellar 
wind. The parameters of the models corresponding to the 
present state of Star 32 are given in Table 1. The variations 
of the stellar wind mass flux with the evolutionary time, for 
different evolutionary sequences, arc shown in Fig. 2. 

The summary of the results can be brief. Equation (5) 
does not seem to correspond to the true history of the stel- 
lar wind strength in LMC X-l/Star 32 binary system (see 
the discussion in Section 3.2). The formulae (6) and (7) (al- 
though obtained ad hoc) seem to be much more realistic, 
and I believe that they give a crude but reasonable descrip- 
tion of the evolution of the stellar wind mass flux from Star 



32. The differences between the results obtained with for- 
mula (6) and formula (7) are not significant (see Table 1 
and Figs Ic and Id), which makes our conclusions more ro- 
bust. To conclude: the present evolutionary mass of Star 32 
must be in the range 31 to 35.5 M©. This result is in very 
good agreement with the new orbital solution of Orosz et 
al. (2009), who give a range for the mass of the star of 28.3 
to 35.3 Mq. According to our evolutionary sequences, the 
initial mass of Star 32 had to be in the range 35 to 40 Mq 
and the present age of the system is 3.7 to 4.0 million Myr. 

In addition, note that if one disregards the estimate of 
M, as a not very reliable constraint, then one is left with 
sccjuciices 1 and 2. The present mass of Star 32 is then in 
the range 32 to 36 Mq, its initial mass is in the range 33.5 
to 37.5 Mq, and the present age of the system is in the 
range 3.9 to 4.2 Myr. As can bo seen, the results are not 
very different. This confirms that our conclusions are not 
sensitive to the uncertainties of the mass-loss treatment. 



4 CONCLUSIONS 

(i) The calculations modelling the evolution of Star 32, un- 
der different assumptions about the shape of the for- 
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time [mln yr\ 

Figure 2. The stellar wind mass loss-rate as a function of evo- 
lutionary time. Ea<;h sequence is labelled with the number corre- 
sponding to that in Table 1. The line no. 9 describes (approxi- 
mately) the rate observed now for Star 32. 

mula describing the stellar wind mass-loss rate, provide 
robust limits on the present mass of the star. For the 
reasonable versions of this formula, this mass should be 
in the range 31 to 35.5 M0. 

(ii) The the initial mass of Star 32 had to be in the range 
35 to 40 Mq. 

(iii) The present age of the system (counted from its birth 
or from the resetting of the evolutionary clock after the 
possible, but unlikely, mass transfer) is in the range 3.7 
to 4.0 Myr. 

(iv) The evolutionary estimate of the mass of Star 32 re- 
mains in very good agreement with the estimate based 
on the orbital solution of Orosz et al. (2009): 28.3 to 
35.3 M0. 
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